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ABSTRACT 

High quality archival spectra of interstellar absorption from C I toward 9 
stars, taken with the Goddard High Resolution Spectrograph on the Hubble Space 
Telescope, were analyzed. Our sample was supplemented by two sight lines, 23 Ori 
and (3^ Sco, for which the C I measurements of Federman, Welty, & Cardelli 
were used. Directions with known CH+ absorption, but only upper limits on 
absorption from C2 and CN, were considered for our study. This restriction 
allows us to focus on regions where CH+ chemistry dominates the production of 
carbon-bearing molecules. 

Profile synthesis of several multiplets yielded column densities and Doppler 
parameters for the C I fine structure levels. Equilibrium excitation analyses, us- 
ing the measured column densities as well as the temperature from H2 excitation, 
led to values for gas density. These densities, in conjunction with measurements 
of CH, CH+, C2, and CN column densities, provided estimates for the amount 
of CH associated with CH+ production, which in turn set up constraints on the 
present theories for CH+ formation in this environment. We found for our sample 
of interstellar clouds that on average, 30-40% of the CH originates from CH+ 
chemistry, and in some cases it can be as high as 90%. 

A simple chemical model for gas containing non-equilibrium production of 
CH+ was developed for the purpose of predicting column densities for CH, CO, 
HCO"*", CH^, and CH.^ generated from large abundances of CH+. Again, our re- 
sults suggest that non-thermal chemistry is necessary to account for the observed 
abundance of CH and probably that of CO in these clouds. 
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1. Introduction 

Carbon is rarely found in the form of C I in the general interstellar medium (ISM). 
In the hot, ionized part of the ISM (like the interior of superb ubbles), it is highly ionized, 
while in the warm neutral regions it is mainly in the form of C II. As the cloud gets colder 
the recombination of C II and the transition from C I to CO appears almost simultaneously 
(HoUenbach & Tielens 1999), leaving little room for C I as the dominant form of carbon. 
The diffuse clouds where this double transition usually occurs are very important for un- 
derstanding the nature of atomic and molecular chemistry in the ISM. It is here where the 
important nonthermal effects, dissipation of turbulence (Duley et al. 1992; Falgarone, Pineau 
des Forets, & Roueff 1995; Hogerheijde et al. 1995; Joulain et al. 1998; Pety & Falgarone 
2000) and MHD waves (e.g., Draine & Katz 1986; Pineau des Forets et al. 1986; Gredel, van 
Dishoeck, & Black 1993; Charnley & Butner 1995; Federman et al. 1996a), can seriously 
alter the chemical composition of the clouds. For example, it is not possible to understand 
the observed abundance of CH+ without taking into account these effects. 

Our first objective is to derive the gas density and the kinetic temperature for such 
diffuse clouds. We infer this information from the observed populations among fine structure 
levels of the electronic ground state of C I. We searched the Hubble Space Telescope (HST) 
archive for spectra of bright stars known to show weak or moderate molecular absorption to 
minimize the effects of denser molecular cores. In all, we found nine suitable target stars; 
their properties from the Simbad database are listed in Table 1. The reddenings are from 
Savage et al. (1977). We included 23 Ori and Sco (not listed in Table 1) in our sample 
with C I column densities from Federman et al. (1997). 

The target stars are located in Taurus, Perseus, Orion, Scorpius, and Ophiuchus, with 
the exception of HD 112244. The general environment toward these stars is well characterized 
through numerous studies, f] Tau, for example, is one of the brightest members of the Pleiades 
(White 1984). Krelowski, Megier, & Strobel (1996) found two absorbing components toward 
the Perseus 0B2 association, but simple molecules (CH, CN) could be detected only in one 
of them. UV observations of Martin & York (1982) and Vidal-Madjar et al. (1982) detected 
several absorption components toward the Perseus association. With the exception of one 
cold component, these were too weak and too hot to be observable in neutral carbon. The 
temperatures deduced by Trapero et al. (1992) from interstellar K I absorption toward stars 
in Perseus are in general agreement with the results of this study. 

de Geus (1992) described the immediate vicinity of the Scorpio-Centaurus region in 
detail. He found that the ISM has been shocked and disturbed by the strong winds of O 
and B stars and past Type II supernovae from the numerous OB associations located here. 
A CO survey of the Ophiuchus region by de Geus, Bronfman, & Thaddeus (1990) found 
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5 main and several weaker cloud components. It is likely, however, that some of them are 
dense molecular cores not relevant for this study. High resolution spectroscopy of Na I D and 
K I A7699 lines by Kemp, Bates, & Lyons (1993) revealed only 3 components, two stronger 
and one weak, in general agreement with our findings of two components. The more recent 
atomic surveys of Welty, Hobbs, & Kulkarni (1994) and Welty & Hobbs (2000) also indicate 
that our C I measurements probe the strongest components toward our sample of stars. 
Copernicus observations of a Scorpii (Allen, Snow, & Jenkins 1990) revealed 4 components 
based primarily on H2 lines. 

The main objective of this study is to obtain information on the importance of non- 
equilibrium CH"*" chemistry for production of CH and CO in diffuse clouds. These molecules 
can be produced in very different environments. In the denser portions of a cloud they 
are involved in a complex thermal chemistry, while in warmer and more tenuous portions 
they are produced via CH'*'. Since any line of sight can involve several clouds and cloud 
sections, the inferred quantities and parameters are likely the result of several environments. 
If one can set up reasonable estimates of the contributions from the different environments, 
further constraints on the physical processes underlying nonthermal CH'*' chemistry can be 
established. 

We follow two approaches in our study. On the one hand, we estimate the amount of 
CH produced in denser regions (with equilibrium chemistry) and establish lower and upper 
limits on CH produced by the molecular chemistry. The difference between our limits and the 
observations should come from non-equilibrium chemistry, and this has to be accommodated 
by any model attempting to explain CH+ chemistry. On the other hand, we calculate CH 
and CO abundances arising from CH'*' in diffuse clouds using the model of Federman et al. 
(1996a) to have a rough estimate of the amount coming from non-equilibrium chemistry. 
As an improvement on the models in Federman et al. (1996a), we use increased effective 
temperatures not only for ions, but also for neutrals directly produced by ion-ion or ion- 
neutral reactions. The effects on CH were noted by Gredel et al. (1993) and described in 
detail by Bucher & Ghnski (1999). 

Details of the data reduction and measurements are discussed in Section 2. The gas 
densities and temperatures were inferred from relative populations obtained by fitting syn- 
thetic spectra to the observations. We discuss this analysis in Section 3 and describe how the 
gas densities and temperatures were inferred from C I excitation in Section 4. We present 
the calculation of molecular abundances in Sections 5 and 6 for the equilibrium and non- 
equilibrium chemistry, respectively. We discuss our result briefly in Section 7 and summarize 
them in Section 8. 
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2. Measurements and Data Reduction 

We retrieved observations for 9 new lines of sight from the HST archive; most were 
acquired at medium resolution (MR; R=-^=20,000) and a handful were taken at high res- 
olution (HR; R= 100,000). These spectra covered several C I multiplets and forbidden lines 
between 1150 A and 1700 A. The datasets that were retrieved from the HST archive and 
used in our study are presented in Table 2. 

We used the standard IRAF STSDAS package to cross-correlate and coadd the subex- 
posures of the archival spectra. The most serious concern was an error in the background 
correction on some of our preprocessed data. We detected this error on all of our HR spectra 
containing the A1260 multiplet of C I. The flux level of a strongly saturated hue, such as 
those of S II and Si II, should be zero at line center; however, in some cases we measured 
negative or small positive values, clearly indicating an erroneous background correction. We 
could then apply an additional correction to compensate for the effect. We never encoim- 
tercd this problem on MR spectra although we had saturated lines in some cases, probably 
because scattered light was more easy to quantify in grating spectra. 

After calibration and data reduction, the components in the neutral carbon lines were 
identified and the continuum of the background star was removed (continuum placement) by 
fitting it to a low order polynomial. The continuum placement turned out to be difficult for 
some multiplets, due to blending of several lines or simply because the lines were saturated. 
For example, in the vicinity of multiplet A1260, Si II and Fe II lines seriously affected the 
continuum placement. 

Following line identification and continuum placement, the Doppler shift (relative to the 
laboratory wavelength) and equivalent width {W\) of neutral carbon lines were measured. 
In some cases, only upper hmits could be extracted from the average Full Width at Half 
Maximum (FWHM) of the instrumental function and the root mean square (rms) average of 
the noise in the stellar continuum (1^;^ < 2 x FWHM x rms). The measured values of Wx 
are listed in Tables 11-19 in the Appendix for comparison with results of other investigations. 

The HST data for e Per and uj^ Sco were supplemented with spectra for additional 
lines acquired with the Copernicus satellite. These coadded spectra were obtained from the 
Multiwavelength Archive at the Space Telescope Science Institute (MAST). All spectra for 
each of the two directions were analyzed together. 
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3. CI Column Densities and Doppler Pcirameters 

The column densities (A^) and the Doppler parameters (6-values) that describe the 
absorbing components along each line of sight were determined by fitting synthetic profiles 
to the observed spectra. The synthetic fluxes were compared to the observations and the set 
of initial parameters was adjusted by a minimization until a satisfactory fit to all observed 
spectra for a sight line could be achieved. Our approach was a modified second order gradient 
method, taken from Bevington & Robinson (1992). The program was tested by reproducing 
the results of curves of growth. We estimated the errors for the inferred parameters to be ~ 
10%, based on the goodness of the continuum placement and the noise in the spectra. The 
effects of other potential sources of error were negligible. For oscillator strengths, we used 
those from the compilation of Morton (1991), updated with our more recent astronomical 
determinations (Zsargo, Federman, & Cardelli 1997; Federman & Zsargo 2001). 

Throughout the synthesis we used a Gaussian instrumental function with a FWHM of 
^ 3 km s^^ and ?a 15 km s^^ for HR and MR spectra, respectively. The exact values of 
FWHM were calculated for each spectrum from the wavelength dependence in the resolution 
of Echelle-A and G160M (described in the instrumental handbook for the GHRS). 

The results of the profile syntheses are shown in Table 3. To demonstrate the quality of 
fits that could be achieved by this procedure, we display observed and synthesized spectra 
for 1 Sco in Figure 1. The differences between the observed and model spectra in the Figure 
are well within the continuum noise for all C I multiplets. The HR spectra of A Ori, 1 Sco, 
5 Sco, and a Sco revealed that there are at least two components in C I absorption along 
these sight fines. The two components toward 1 Sco, for example, are clearly visible in the 
HR spectra of C I A1328.83 in Figure 1. We could not discern multiple components toward 
most sight lines for which only MR spectra were available. The only exception is HD 112244 
where the large separation (^ 24 km s~^) between the components was easily recognizable. 
The main absorption component (Component A) for the lines of sight in Table 3 are very 
well defined. The uncertainties in their columns and 6-values are dominated by the errors 
in the continuum placement and by the signal-to-noise ratio of the spectra. The parameter 
values for the secondary components (Component B), on the other hand, are not defined 
very well. Often, we had to keep some parameters fixed during the synthesis in order to 
achieve reasonable values for the column densities, or else, only upper limits of N could be 
provided. 

Comparison between our line parameters and those extracted from high and ultra-high 
resolution measurements of atoms and molecules shows good correspondence. There are four 
sight fines with one C I component, rj Tau, e Per, tt Sco, and lo^ Sco. The sight fines toward 
rj Tau and e Per indicate but one component in K I (Welty & Hobbs 2000) with 6-values 



-6- 



of about 1 km s~^. Presumably, the larger 6- values in C I arise because this atom probes 
more widely distributed material. For n Sco and Sco, Welty et al. (1994) and Welty 
& Hobbs (2000) see multiple components in Na I D and K I, respectively. Here, our lower 
resolution measurements could not distinguish absorption from the stronger, closely spaced 
components. The 6-values for the J = lines reflect this. The narrower 6-values for lines 
from J = 1 and 2 are likely the result of their sampling only the strongest component, a 
smaller portion of the hue of sight, or both. 

A similar pattern emerges for the sight lines with two C I components. The ^ 12 km s^^ 
separation seen toward A Ori comes from the complexes discerned in Na I D measurements 
(Hobbs 1974; Welty & Hobbs 2000). Toward HD 112244, the two components about 24 km 
s~^ apart are also present in CH spectra (Banks, Federman, & Lambert 1984). For 1 Sco, 
the correspondence is seen in Na I D complexes (Welty et al. 1994; Welty & Hobbs 2000), 
separated by about 7 km s~^. The same apphes to the two C I components (A^; Ki 5 km s~^) 
toward 6 Sco, which are present in Na I D and K I (Welty & Hobbs 2000) and CH"*" (Crane, 
Lambert, & Sheffer 1995). Finally, the 9 km s^^ separation in components toward a Sco is 
seen in Na I D and K I (e.g., Welty & Hobbs 2000). 



4. Excitation Analysis 

We used the column densities found from the profile syntheses to infer gas densities and 
temperatures for the gas containing C I. To do so we assumed steady-state detailed balance 
among the fine structure levels J = 0, 1, and 2 of neutral carbon. Two independent number 
density ratios can be estabhshed. These are 

nj=2 _ -Ro,2 (^1,0 + -Rl,0 + -^1,2) + -Rl,2-Ro,l 

~ (A2,l + R2,l + R2,o) (Ai^O + Rl,0 + Rl,2) - Rl,2 (A2,l + i?2,l) 

and 



nj=i ^ Ro,i (^2,1 + R2,l + R2,o) + Ro,2 {A2,l + R2,l) 

nj=0 (^2,1 + R2,l + i?2,o) (^1,0 + Rl,0 + Rl,2) - Rl,2 (^2,1 + i?2,l) ' 

where the terms Rij are defined by 



(2) 



k 

The quantities are the coUisional excitation or de-excitation rate coefficients through 
collision with species k and Gij are the ultraviolet radiative pumping rates between fine 
structure levels i and j. In eqns. 1 and 2, the quantity A^^i is the far infrared spontaneous 
decay rate from level i to level 1. Variables Uk are the number density of species k. For our 
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lines of sight, collisions with H, He, and H2 were considered, since they are the most abundant 
species in diffuse interstellar clouds. The C I coUisional de-excitation rate coefficients are 
listed in Table 4. The corresponding excitation rates can be calculated by the formula 



2z + 1 
'2j + l' 

where kTij and T are the energy between states i and j and the gas temperature, respectively. 
The rates for C I collisions with H2 and H were taken from Schroder et al. (1991) and from 
Launay & Roueff (1977), respectively. The coUisional rates for He were from Lavendy, Robbe, 
& Roueff (1991) and Staemmler & Flower (1991). Radiative pumping was neglected in our 
analysis because it has a negligible effect on the populations, unless the ultraviolet field is 
at least 10 times the strength of the average interstellar field. 

If the clouds have close to uniform density and chemical composition, then 

Uk ~ , 5 



nj=i _ Nj=i nj=2 _ Nj=2 

^ (XTld ^ , 

nj=o Nj=o nj=o Nj=o 

where the A'"-s are the column density of species k, hydrogen, H2, hehum, or a C I fine 
structure level, and n is the gas density. We calculated the fractional abundances for H 
and H2 using column densities from Savage et al. (1977) and we assumed that 10% of the 
material is He for all lines of sight. 

After substituting eqn. 5 in eqns. 1 and 2, one has two equations for the unknown 
temperature and gas density. Unfortunately, the temperature and density solutions for 
eqns. 1 and 2 are not unique. One can find many density-temperature pairs to satisfy these 
equations for a specific ^^^^^ or ^^^^ ratio. As a further complication, the uncertainties in 
the C I column densities introduce uncertainties in the C I level ratios; therefore, a density 
range rather than a given value satisfies eqn. 1 or 2 for a temperature. As a consequence 
of this ambiguity, numerous density-temperature pairs simultaneously satisfy eqns. 1 and 2. 
One needs independent density or temperature estimates to further confine the solutions. 

The density-temperature solutions of eqns. 1 and 2 are displayed in Figure 2 for the main 
components of each sight fine. The C I excitation analysis was supplemented by the temper- 
ature estimates of Savage et al. (1977) which were determined by H2 excitation analysis for 
the sight fines. Since the H2 temperatures of Savage et al. (1977) are valid only for the main 
components, we did not perform the C I excitation analysis for the secondary components. 
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The sole exception is HD 112244 where the two components are fairly similar; therefore, the 
temperature estimate of Savage et al. (1977) most likely applies to both components. In 
Figure 2, each C I ratio resulted in two density-temperature curves corresponding to the 
maximum and minimum ratios that were allowed by the uncertainties in the C I column 
densities. Any density-temperature point between the solid or the dashed curves in Figure 2 
is consistent with the observed C I J=0 and 1 or the J=0 and 2 column densities. Then, 
the region that is between all four curves contains all density-temperature pairs that are the 
simultaneous solutions to eqns. 1 and 2. However, Figure 2 also shows that the temperature 
values are not well defined for most sight lines. The allowed temperatures range between 30 
and 200 K, and can be confined only by the estimates of Savage et al. (1977). For simplicity, 
we accepted as solutions all the density-temperature pairs in Figure 2 that lie within the 
thick rectangles. These regions provide the smallest independent density and temperature 
ranges that satisfy both eqns. 1 and 2 for the observed C I ratios and consistent with the H2 
excitation temperatures. The derived density and temperature ranges are listed in Table 5. 
In the case of rj Tau, we chose to list only an upper limit for the density, since the lower 
limit is essentially zero. 



5. CH Production in Moderately Dense Regions 

To estimate the amount of CH produced in the denser regions in our sight lines, we 
followed the method described in detail by Federman et al. (1997). In their work, they only 
used the most important processes to find the analytic expression 

"^^^"^ - G(CH) 

that describes the CH fractional abundance. In eqn. 6 and all the following equations, ki 
denotes the rate coefficient for reaction i and the fractional abundance of a species y is x{'y)— 
n(y)/(n(H) -|- 2n(H2)). The terms n, n(y), and G(y) are the gas density, number density of 
species y, and the photodissociation rate for species y, respectively. The photodissociation 
rate includes grain attenuation of ultraviolet radiation which is approximated by tuv = '^Ay 
(Federman et al. 1994). For the sight lines in Scorpius, which have larger ratios of total to 
selective extinction, a value of lAAy was adopted. One can invert eqn. 6 and solve for n 
using the observed CH abundance for a comparison between the densities derived from C I 
excitation. We can also constrain our results for CH, using available upper limits for N{C2) 
and A^(CN). In equilibrium, the fractional abundances are given by the following expressions 
(Federman et al. 1997): 



(7) 
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and 



-e^^) - (air) ( 



kix{C+) 




The rate coefficients and photodissociation rates used in this analysis are given in Table 6. 
A detailed description of these coefficients is given by Federman et al. (1994) with updates 
in Pan, Federman, & Welty (2001). Table 7 lists the necessary fractional abundances and 
Tuv calculated from E{B — V) values of Savage et al. (1977). The data for H2 and the 
carbon-bearing molecules, CH, C2, and CN, are from Savage et al. (1977) and Federman 
et al. (1994), respectively. The CO abundances are from recent analyses of Copernicus 
spectra (Crenny & Federman 2001, private communication). We used the CH+ abundances 
of Federman (1982), Lambert & Banks (1986), and Price et al. (2001). The data for x{C+) 
are from Sofia et al. (1997), x{0) from Meyer, Jura, & Cardelli (1998), and x(N) from 
Lugger et al. (1978) and York et al. (1983). 

We estimated the gas densities necessary to reproduce the observed CH columns in an 
environment in thermal equilibrium by inverting eqn. 6 and using the data in Tables 6 and 7. 
Such environments are very different from those required to produce the observed amount of 
CH+. The calculated upper limits on densities, as expected, are much larger than those found 
from C I excitation (see Section 5). The sole exception is a Sco, where the two densities are 
marginally consistent with each other. Since we chose lines of sight with no detectable C2 and 
CN absorption, densities of several thousand cm~^ are not realistic, and the values inferred 
from C I excitation are more appropriate. The clear inconsistency resulting from the two 
density calculations is a strong indication that non-equilibrium CH chemistry is important 
toward these stars. This contrasts with the results for the molecule-rich gas toward ( Oph 
(Lambert et al. 1994) where the above calculation and the observations gave similar results. 

Using a;(C2) or a;(CN) as well as Ntotai from Table 7, we can estimate the column densi- 
ties of CH associated with the equilibrium chemistry from eqns. 7 or 8. Then, a lower limit 
on the amount of CH associated with the non-equilibrium CH+ chemistry can be imposed 
by comparing the calculated amounts to the observed column densities. On the other hand, 
if we use the densities derived from C I excitation in eqn. 6, an upper limit on the amount of 
CH produced via CH"*" can be obtained. In Figure 3 we display the results, normalized to the 
observed column of CH+. The triangles and downward-pointing arrows represent the mea- 
sured values and upper hmits of the A'"(CH)/A'"(CH+) ratios, respectively. Ranges indicate 
the normalized CH columns that are likely to be associated with the non-equilibrium CH"*" 
chemistry. Again, the results in Figure 3 strongly suggest that non-equilibrium chemistry is 
indeed very important for the production of CH in molecule-poor sight lines, and in some 
cases, like that of 23 Ori, it is the dominant source for this molecular species. 
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A larger contribution from the non-equilibrium CH"*" chemistry would be necessary if 
CH is effectively destroyed by reactions with atomic hydrogen (reaction 30 in Table 8) in the 
moderately dense regions along our sight lines. This reaction is thought to be suppressed 
by a reaction barrier (e.g., Mitchell & Deveau 1983); however, recent quantal calculations 
(Harding, Guadagnini, & Schatz 1993; van Harrevelt, van Hemert, & Schatz 2002) suggest 
that this may not be the case. Unfortunately, the rate coefficient for the reaction is not well 
defined for the gas temperatures relevant to the diffuse ISM, and the value of a;(H) in the 
moderately dense regions is not known — only line of sight averages are available. Therefore, 
the importance of the reaction, CH + H — > C + H2, could not be judged properly and the 
reaction was not included in eqns. 6 — 8. If significant amounts of CH are destroyed through 
this reaction, our conclusion about the importance of non-equilibrium chemistry would be 
strengthened. 



6. Nonthermal Chemistry in Diffuse Clouds 

We used the chemical model of Federman et al. (1996a) to examine the effects of non- 
thermal motions on diffuse cloud chemistry. The origin of the nonthermal motion is assumed 
to be the result of the propagation of Alfven waves, constantly entering the cloud from the 
intercloud medium. The dissipation of these waves is concentrated primarily in a diffuse 
cloud boundary layer, the transition zone between the cloud and the intercloud medium. 
The molecular abundances resulting from the Alfvenic driven nonthermal chemistry are cal- 
culated for this boundary layer, without considering any time dependence or attempting to 
describe any changes in physical or chemical conditions throughout the region. Instead, we 
impose a set of conditions for the clouds toward our sample of stars that are inferred from 
previous observations or are from the excitation analysis of neutral carbon. 

In the chemical model we use a set of 16 atomic and molecular species closely related 
to CH+ chemistry. These are H, He, C, O, He+, C+, H2, CO, CH, OH, H2O, C0+, CH+, 
CH^, CH^, and HCO+. The set of molecules and chemical reactions was chosen so that 
it constitutes a closed system with reasonable assumptions on the physical conditions in 
diffuse clouds. The chemical reactions involved in this model are given in Table 8; the rate 
coefficients were taken from Federman et al. (1997) and Drainc & Katz (1986), updated 
with information from van Dishoeck & Black (1988), Federman & Huntress (1989), Mitchell 
(1990), Anicich (1993), Millar, Farquhar, & Willacy (1997), and Pan ct al. (2001). For 
example, the CO photodissociation rates come from the tabulation in van Dishoeck & Black 
(1988) using the observed column densities of CO and H2 and tuv- The photodissociation 
rates for the other molecules included grain attenuation through the tuv values listed in 
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Table 7. Abundances of CH, CO, C0+, HCO+, He+, CH+, and CH+ were calculated by 
solving a system of nonlinear rate equations. Wc assumed that for all lines of sight the 
He, OH, and H2O fractional abundances are 0.1, 1.6 x 10~^°, and 3.2 x 10^^^, respectively, 
based on the theoretical results of Federman et al. (1996b). The abundances of the remaining 
elements and molecular species, listed in Table 7, were taken from observations (see Section 5) 
and were kept fixed throughout the analysis. Fractional abundances of neutral carbon are 
from the present analysis. The electron densities were derived from the 0+ abundances. 

Wc define an effective temperature for each ion-neutral (and in some cases for each 
neutral-neutral) reaction in order to introduce the effects of nonthermal motions. Federman 
et al. (1996a) defined T^ff as 

h^ff = hr + ^i^vl^ude ■ (9) 

where k and // are the Boltzmann constant and the reduced mass of the system, respectively. 
The term Vpartide is the turbulent velocity of either the ionic (vion) or the neutral {vneutrai) 
species and defined by the root mean square Alfven wave amplitude or zero, depending on 
the physical environment wc attempted to model. 

First, wc tested our code by attempting to reproduce the results of Federman et al. 
(1996a). They imposed a set of conditions on their model clouds that are believed to be 
representative of the region in which the effects of nonthermal chemistry are significant. 
Therefore, we modeled their "hypothetical" cloud with fractional abundances of 0.2, 0.076, 
2 X 10-^ 7.94 X 10-^ and 3.02 x 10"^ for H2, He, C+, N, and O, respectively The 
visual extinction was set to 0.3 mag, representative of a predominantly atomic diffuse cloud, 
corresponding to a value of 0.6 for tuv- Since we did not calculate CH+ abundances, as 
opposed to Federman et al. (1996a), we used their result for CH"'" to calculate the amounts of 
CO"*", HCO"*", and CH produced by the nonthermal chemistry. The gas density, temperature, 
and Vion combinations used for the modeling correspond to their models 1, 2 and 3 — (100 
cm-3, 100 K, 3 km s'^), (50 cm~^ 100 K, 3 km s'^), and (200 cm^^ 100 K, 3 km s'^), 
respectively. In all cases the calculated CO"*", HCO"*", and CH abundances were within 
30-40% of the values listed in Table 1 of Federman et al. (1996a). 

In the models of Federman et al. (1996a), the nonthermal effects were considered only for 
charged species, reffecting the distinction that ions are tied to the magnetic field. However, if 
one considers the conservation of momentum during a chemical reaction, then it is plausible 
that the product species (ion or neutral) could carry some of the momentum from the parent 
molecules (Gredel et al. 1993; Bucher & Ghnski 1999). This would mean that even neutral 
species, directly produced by reactions involving ionic species, can have more turbulent 
velocity than the average value for neutrals (not related directly to the ionic species). One 
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could establish a complex model, where the turbulent velocity progressively decreases from 
the value representative for ions to that of the average for neutrals, following the chain of 
chemical reactions from ionic to neutral species. Such an approach is beyond the scope of 
this study. Instead, we simply explored the phenomena by setting the turbulent velocity for 
neutral molecules either to the ionic turbulent velocity or zero. 

In our approach, we modeled three different situations, one with no turbulence at all 
(e.g. equilibrium chemistry), one with Vion 7^ but Vneutrai = 0, and one with nonzero 
Vion and Vneutrai = Vion- Wc uscd thc equilibrium CH+ abundance of 3.3x10^^^ calculated by 
Fcdcrman ct al. (1996a) whenever we attempted to model equilibrium chemistry without any 
turbulence. Since the observed CH+ abundance is primarily the product of non-equilibrium 
chemistry, we would introduce the effects of this chemistry into our "equilibrium" run {vion 
— and Vneutrai — 0) by using the observed amount of CH+. However, we used the observed 
amount of CH+ in the cases when the turbulent velocities were set to nonzero. The values 
for turbulent velocity were chosen to be the amplitude of Alfvenic waves. We adopted 
a standard value of 3 km s^^ for both Vion and Vneutrai (whenever they differ from zero), 
following Federman et al. (1996a), based on the typical width of CH^ absorption lines. We 
used the temperature and density from the C I excitation analysis (see Table 5) as T and n. 

The results of our calculations arc shown in Tables 9 and 10 for CH and CO. It is 
immediately obvious that the observed abundances cannot be reproduced by equilibrium 
chemistry alone; there arc several orders of magnitude difference between the observed and 
calculated abundances for all lines of sight for both CH and CO. Our results indicate that a 
nonthermal chemical model for ionic species only produces too much CH. When we also take 
nonthermal effects into account for neutral-neutral reactions, the observed CH abundances 
are reproduced reasonably well. The lone exception is the unusual line of sight toward 23 
Ori, which has detectable amounts of CH and CH+ even though the H2 column is about 10^^ 
cm~^. When Tgff increases, reaction 17 (see Table 8) becomes more effective in destroying 
CH. This reaction also increases the amount of CO between models 2 and 3. Unlike the 
situation for CH, we did not find reasonable agreement for CO in general, indicating the 
limits of our approximations. Thc closest correspondence occurs for the denser sight lines, 
those in Scorpius, suggesting that pockets of higher density control the abundance of CO in 
diffuse gas. Since a theoretical determination of the rate coefficient for 0+ + OH — ^ C0+ + 
H (Dubernet, Gargaud, & McCarroU 1992) suggests a value much larger than that used here, 
we investigated its affect on CO as well. While the equilibrium abundance for CO increases 
nearly 5 fold, the non-equilibrium results are changed only 5%. Clearly, the production of 
CO^ and HCO"*", the precursors of CO, is dominated by the presence of enhanced amounts 
of CH"*" brought about by non-equilibrium chemistry in thc sight lines considered here. It is 
likely that we have to improve our model in other ways as well. For instance, a more precise 
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description for Vneutrai is needed as is a more detailed analysis of H2 and CO photodissociation 
that includes self- and mutual-shielding. 

7. Discussion 

Observations and models have suggested long ago that the ISM is not in an equihbrium 
state. The underlying reason is that interstellar clouds cannot be treated as a closed systems; 
interactions with often violent surroundings are inevitable. Shock fronts from supernovae 
and cloud collisions are probably commonplace. These non-equilibrium processes affect the 
chemistry of most molecular species, including that of CH"*". In fact, almost everyone agrees 
that these non-equilibrium processes hold the key for explaining the elusive CH^ abundance 
problem in the ISM. The exact mechanisms and place for CH+ production, however, is still 
the subject of intense debate. MHD shock models (e.g., Flower & Pineau des Forets 1998) 
or those that involve intermittent dissipation of interstellar turbulence (e.g., Joulain et al. 
1998) suggest that CH+ is formed in warm material at temperatures T >100 K. Excitation 
analysis of interstellar C2 (Gredel 1999), on the other hand, places the production of CH+ 
into the colder gas at temperatures of 50-100 K. Our results are more consistent with this 
later suggestion and support the claim that non-equilibrium chemistry is very important 
(or sometimes even dominant) for CII+, CH, and CO in diffuse sight lines. We concur 
with Credell (1999) that about 30% of the observed CH is associated with CII+ synthesis. 
The associated molecular ions, CO"*" and HCO+, are likely to be affected as well, with 
respective abundances reaching ~ 6 x 10~^^ and about 5 x 10~^^. (We note in passing 
that such enhancements in the HCO"^ abundance still yield column densities at least an 
order of magnitude below those inferred from absorption of the millimeter- wave continuum 
from extragalactic sources [Lucas & Liszt 1996].). This conclusion highhghts the fact that 
any model attempting to account for CH"*" abundances must address the "side-effects" on 
other species and be able to describe their abundances. Therefore, indirectly, we imposed 
constraints on the suggested theories thought to be responsible for CH"*" production in this 
environment. 

There were cases where our simplified models failed, indicating some shortcomings in 
our approach. It remains to be seen whether the source of these discrepancies is related to the 
assumed physical environment for particular lines of sight (gradients in physical conditions) 
or to effects ignored in our approach, such as detailed calculations for II2 and CO photodis- 
sociation. Another possible source of error is the "uniformity" of our approximations. The 
differences in 6-values and Doppler shifts associated with the C I fine structure levels for 
a given sight hue indicate that the levels might not fill the same volume. The excitation 
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may occur in clumps. A more careful excitation analysis of neutral carbon is necessary to 
infer average density and temperature for diffuse clouds, but the means to perform such an 
analysis are beyond the capabilities of current theory and observation. Moreover, we either 
applied a given turbulence for neutrals or we did not. In reality, the turbulence can change 
gradually from the values representative for ions to that of neutrals as one follows the chem- 
ical network from an ionic parent molecule (like CH"*" in our case) to products (Bucher & 
Ghnski 1999). Further work should explore the importance of these various effects. 

8. Summciry 

We surveyed the physical and chemical conditions in the diffuse ISM with the aim of 
assessing the importance of non- equilibrium processes. Wc retrieved a mixture of high and 
medium resolution spectra from the HST archive. The selection criteria for the lines of sight 
was motivated by the idea to minimize the contribution from the denser diffuse cloud cores, 
which are satisfactorily described by equilibrium chemistry. We studied transitions involving 
the fine structure levels of neutral carbon to infer the physical parameters of these clouds 
through analysis of excitation. Our conclusions can be summarized as follows: 

We could identify multiple cloud components along the sight fines for which high- 
resolution spectra were available and toward HD 112244 where the separation between the 
components was larger than the resolution of the MR spectra. We derived C I column 
densities and 6-values for the components or provided upper limits for non-detections. 

We performed C I excitation analyses for all main components to infer n and T. The gas 
densities and kinetic temperatures were comparable to the representative values for diffuse 
clouds. 

We calculated the abundances produced by molecular chemistry in the moderately dense 
regions along our sight lines. By removing the contribution from the molecular chemistry, 
we quantified the amounts of CH from nonthermal chemistry. Our analysis suggested that 
the amount of CH associated with CH+ is at least 30-40% of the amount of CH in general, 
but in some cases it can reach as high as 90% (e.g., 23 Ori). 

Density estimates based on the observed upper fimits for C2 and CN, along with CH 
columns, and on the assumption that the equilibrium chemistry acts alone were (much) larger 
than the densities from C I excitation. The discrepancy between these density estimates 
clearly indicates that the assumption of equilibrium chemistry as the sole contributor is 
inappropriate. 
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We modeled the non-equilibrium chemistry related to CH+ and improved upon the 
approximation of Federman et al. (1996a) to account for nonthermal motions, the driving 
force behind non-equilibrium chemistry. We defined effective temperatures for both ions and 
neutrals which involved thermal and nonthermal motions. We explored the effect of these 
motions with different combinations of neutral and ion turbulent velocities. We found that 
the observed amounts of molecular species cannot be reproduced by equilibrium chemistry. 
Models with non-zero turbulence for ions only also fail; turbulence for both ionic and neutral 
species was required to match the observations satisfactorily. 

In closing, non-equilibrium chemistry is needed to explain the physical and chemical 
conditions in the diffuse gas toward our sample of stars and substantial amounts of CH are 
produced by nonthermal chemistry in the diffuse ISM along these sight lines. 

A. Appendix 

Here we present tables of equivalent widths for C I hues seen in our spectra (Tables 11- 
19). Comparison is made with other determinations, when available, and with the results of 
our profile syntheses. The agreement among the various measures is quite good, especially 
when one considers the fact that the MR spectra from HST and the Copernicus spectra do 
not always resolve blended lines or velocity components. Specific examples are given in the 
notes to each table. 

We thank Dan Welty for his comments on the Ph.D. Dissertation from the Univ. of 
Toledo on which this paper is based. Comments by the anonymous referee improved several 
points in our paper. This research was supported by STScI grant AR-08352.01-97A. For the 
Copernicus data, we utilized the archive developed by George Sonncborn and available at 
the Multiwavelength Archive at STScI. We made use of the Simbad database, operated at 
Centre de Donees Astronomiques de Strasbourg, Strasbourg, France. 
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Table 1. Data for Target Stars 



HD 


Star 


Galactic Coordinates 


Distance 


Spectral type 


V 


E{B - V) 






in 




(pc) 




(mag) 


(mag) 


23630 


T] Tau 


166.67 


-23.46 


113 


B7IIIe 


2.87 


0.00 


24760 


e Per 


157.35 


-10.09 


165 


B0.5V 


2.90 


0.09 


36861 


A Ori 


195.05 


-12.00 




08III 


3.30 


0.12 


112244 




303.55 


6.03 


578 


09Ibe 


5.32 


0.34 


141637 


1 Sco 


346.10 


21.71 


160 


B1.5Vn 


4.64 


0.20 


143018 


TT Sco 


347.21 


20.23 


141 


BIV 


2.89 


0.08 


143275 


6 Sco 


350.10 


22.49 


123 


B0.2IVC 


2.29 


0.16 


144470 


ui^ Sco 


352.75 


22.77 


130 


BIV 


3.95 


0.22 


147165 


a Sco 


351.31 


17.00 


225 


BlIII 


2.91 


0.38 
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Table 2. HST Observations for Our Sight Lines 



Star 


Dataset 


T] Tau 


zlbw020st,zlbw020tt,zlbw020ut,zlbw020vt,zlbw020wt, 




zlbw020yt,zlbw0210t 


e Per 


z2c20206t,z2c20207t,z2c20208t 


A Ori 


Z3di0204t,z3di0205t,z2c20306t,z2c20307t,z2c20308t, 




z2520306t.z2520307t,z2520308t 


HD 112244 


z0yv0506m,z0yv050gt,z0yv060em 


1 Sco 


Z2zx020bp,z2zx020fp,z0zi0708m,z0zi070bm,z0zi070ct, 




zOziOTOdt ,z0zi070gt ,z0zi070it 


TT Sco 


Z0zi0408t,z0zi0409t,z0zi040at,z0zi040bt,z0zi040ct, 




Z0xz0308t,z0xz0309t,z0xz030at,z0zi040et,z0zi040gt 


5 Sco 


Z3di0504t,z3di0505t,z34r020dt,z34r020et,z2c20506t, 




z2c20507t,z2c20508t 


uj^ Sco 


z2c20606t,z2c20607t,z2c20608t 


a Sco 


Z34r0409t,zl3v0108m,zl3v0109m,zl3v010bm,zl3v010cm, 




zl3v010dni,zl3v010gm 



Table 3. Derived Parameter Values 



Star Component A Component B 

Level Column density fe- value Doppler shift" Column density 6- value 

(cm~^) (km s^^) (km s^^) (cm~^) (km s^^) 






2.4(±0.2) 


X 10^^ 


2.5±0.3 












1 


< 7.7 X 


10" 














2 


< 5.2 X 


10" 














e Per 





















4.4(±0.4) 


X 10^^ 


3.0* 












1 


6.4(±0.6) 


X 1012 


2.1±0.2 












2 


< 8.2 X 


10" 














A Ori 





















1.3(±0.1) 


X 10^^ 


3.0±0.3 


-12.1* 


9.6(±1.0) 


X 


1012 


3.0* 


1 


5.5(±0.6) 


X 10^3 


2.4±0.2 


-12.1* 


1.7(±0.2) 


X 


1012 


3.0* 


2 


9.2(±0.9) 


X 10^2 


2.5* 


-12.1* 


1.4(±0.1) 


X 


1012 


3.0* 


HD 112244 





















2.2(±0.2) 


X 10^^ 


4.5±0.5 


-23.6±0.2 


1.9(±0.2) 


X 


10" 


4.0±0.4 


1 


3.2(±0.3) 


X 10^3 


2.2±0.2 


-24.1±1.7 


3.2(±0.3) 


X 


1013 


2.7±0.3 


2 


4.4(±0.4) 


X 1012 


2.0* 


-23.0* 


6.4(±0.6) 


X 


1012 


2.0* 


TT ScO 





















9.2(±0.9) 


X 1012 


4.4±0.4 












1 


1.4(±0.1) 


X 1012 


2.0* 












2 


< 2.4 X 


10" 














1 Sco 





















2.5(±0.3) 


X 1013 


1.8±0.2 


-7.1* 


2.2(±0.2) 


X 


1013 


3.1±0.3 


1 


3.3(±0.3) 


X 1013 


1.5±0.2 


-7.1* 


4.3(±0.4) 


X 


1012 


2.5±0.3 


2 


1.6(±0.2) 


X 1013 


1.2±0.1 


-7.1* 


4.0(±0.4) 


X 


10" 


2.5* 


6 Sco 





















8.9(±0.9) 


X 1013 


2.3±0.2 


-4.8* 


1.3(±0.1) 


X 


1013 


1.4±0.1 


1 


5.8(±0.6) 


X 1013 


1.2±0.1 


-4.8* 


3.8(±0.4) 


X 


1012 


1.5* 


2 


1.2(±0.1) 


X 1013 


1.2±0.1 




< 6.7 X 


10" 
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Table 3 — Continued 



Star Component A Component B 

Level Column density 5-value Doppler shift" Column density 6-value 

(cm~2) (km s~i) (km s~^) (cm"^) (km s~^) 



u}^ Sco 



















1.6(±0.2) 


X 


1014 


3.0^ 








1 


5.3(±0.5) 


X 


10" 


1.0^ 








2 


1.0(±0.1) 


X 


10" 


LO*" 








a Sco 



















1.1(±0.1) 


X 


10" 


2.7±0.3 


-9.0±0.1 


8.4(±0.8) 


X 10^2 4.0±0.4 


1 


9.9(±1.0) 


X 


10" 


1.8±0.2 




< 2.1 X 


1012 


2 


2.5(±0.3) 


X 


10" 


1.4±0.1 




< 1.2 X 


10" 



^ The quantity "Doppler shift" is a shift relative to Component A. It was measured in 
HR (FWHM w 3 km s~i) spectra and was sometimes adjusted during profile synthesis to 
improve the fit. In the case of HD 112244, the large shift between the two component could 
be measured in MR (FWHM w 15 km s~i) spectra. 

^ The parameter was kept fixed during the fitting procedure. 
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Table 4. CoUisional deexcitation coefficients 



Coeff. Transition from level i to level j Notes 

(cm^ s-i) J=2 ^ J=2 ^ 1 J=l ^ 





7.244 • 


10- 


-11 rpO.069 


1.413 • 


10- 


-10 rpo.ise 


1.349 


• 10 


-10 


'pO.043 


1 


"J 


1.995 • 


10" 


-11 rpO.181 


2.951 • 


10" 


-11 rpO.241 











2,3 


ortho—H^ 


3.311 


• 10- 


-11 rpo.ie 


4.365 


• 10- 


-11 'pO.27 


8.511 • 


10- 


-11 


'j'-O.OS 


4 


para-H2 
^i,3 


1.072 • 


10- 


10 -^-0.04 


1.479 


■ 10- 


-10 rpO.OS 


1.122 • 


10- 


-10 


-p-O.lO 


4 



The quantity T is the gas temperature. 

Notes. - (1) Launay & Roueff 1977; (2) Lavendy et al. 1991; (3) Staemmler & Flower 1991; (4) 
Schroder et al. 1991. 
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Table 5. Gas Densities and Temperatures 



Star 


nc i" 








(cm-3) 


(K) 


(cm-3) 


r] Tau 


<3.5 


59-143 




e Pcr 


11.5-33 


51-188 




A Ori 


68-182 


36-60 




HD 112244 A'^ 


14-28 


67-113 


< 960 


HD 112244 B<^ 


16.5-48 


67-113 


< 960 


1 Sco 


200-800 


53-114 


< 2299 


TT Sco 


13-33 


56-144 


< 2793 


S Sco 


90-265 


42-81 


< 1115 


LU^ Sco 


43-83 


59-94 


< 643 


a Sco 


137-345 


50-87 


< 371 


23 Ori 


7.5-23 


65-150 


< 4004 


Sco 


35-205 


77-103 


< 770 



Density inferred from the C I level populations. 
^ H2 excitation temperatures of Savage et al. 1977. 
" Density calculated by inverting eqn. 6. 

A is the stronger and B is weaker component. 
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Table 6. Chemical Network for the Equihbrium Chemistry 



Reaction Rate coefficient Value Notes 

or Rate 



c+ + H2 ^ cut + ki 

C+ + CH ^ + H k2 

C2 + N ^ CN + C k3 

CH + N ^ CN + H ki 
^ , r CN+ + H 
^"+^^-^1 CH++N 

CH + W ^ C + H G(CH) 

C2 + hi^ ^ 2C G(C2) 

CN + hi/ ^ C + N G(CN) 



1.0 


X 


10" 


cm^ 


1 


5.4 


X 


10- 




2 


1.7 


X 


10- 


V 300 ^ " 


3 


2.0 


X 


10- 


cm^ s-i 

Y 300 ^ " 


4 


5.0 


X 


10- 


10 cm3 s-i 


3 


1.3 


X 


10- 


■9 Q-TUV g-l 


2 


2.0 


X 


IQ- 


-10 g-i"c/v g-1 


5 


1.0 


X 


10- 


-10 Q-ruv g-1 


6 



Notes.- (1) Federman et al. 1997; (2) Pan et al. 2001; (3) Federman et al. 1994; (4) Messing 
et al. 1981; (5) Pouilly et al. 1983; (6) Lavendy et al. 1987. 



Table 7. Fractional Abundances 



a;(Y) 



Star 


TUV 


AT b 

J^total 


C+ 


C 


N 





H 


2 


CH+ 


CO 


CH 


C2 


ON 


r) Tau 


0.00 


2.7(20)^ 


1.4( 


-4) 


9.1( 


-8) 


7.5( 


-5) 


4.0( 


-4) 


2.6( 


-1) 


7.0(-9) 


<8.1(-9) 








e Per 


0.56 


3.2(20) 


1.4( 


-4) 


1.5( 


-7) 


3.2( 


-5) 


3.0( 


-4) 


2.1( 


-1) 


<3.4(-9) 


2.7(-9) 


6.2(-9) 






A Ori 


0.74 


6.3(20) 


1.3( 


-4) 


3.2( 


-7) 


7.9( 


-5) 


3.2( 


-4) 


4.0( 


-2) 


<l.l(-9) 


6.2(-9) 








HD 112244 A 


1.05 


7.5(20) 


1.4( 


-4) 


3.8( 


-7) 


7.5( 


-5) 


4.0( 


-4) 


1.9( 


-1) 


1.0(-8) 


1.8(-8) 


4.0(-9) 




<4.2(-10) 


HD 112244 B 


1.05 


7.5(20) 


1.4( 


-4) 


3.0( 


-7) 


7.5( 


-5) 


4.0( 


-4) 


1.9( 


-1) 


<3.0(-9) 


5.1(-9) 


2.8(-9) 




<4.2(-10) 


1 Sco 


0.87 


1.6(21) 


1.6( 


-4) 


6.4( 


-8) 


1.0( 


-4) 


3.5( 


-4) 


2.0( 


-2) 


<8.8(-10) 


<1.6(-9) 


<l.l(-9) 




<2.4(-10) 


TT Sco 


0.35 


5.6(20) 


1.4( 


-4) 


1.9( 


-8) 


7.5( 


-5) 


3.5( 


-4) 


7.0( 


-2) 


<2.0(-9) 


1.5(-9) 


<8.9(-10) 


<3.0(-9) 




6 Sco 


0.69 


1.4(21) 


1.8( 


-4) 


1.2( 


-7) 


9.8( 


-5) 


4.2( 


-4) 


4.0( 


-2) 


2.4(-9) 


2.2(-9) 


1.6(-9) 




<9.3(-ll) 


ui^ Sco 


0.96 


1.7(21) 


1.4( 


-4) 


1.3( 


-7) 


1.2( 


-4) 


4.3( 


-4) 


1.3( 


-1) 


3.6(-9) 


5.3(-9) 


1.9(-9) 


<1.0(-9) 




a Sco 


1.65 


2.3(21) 


1.4( 


-4) 


1.0( 


-7) 


1.0( 


-4) 


4.3( 


-4) 


5.0( 


-2) 


2.5(-9) 


1.3(-9) 


1.3(-9) 




<5.7(-ll) 


23 Ori 


0.68 


5.5(20) 


2.8( 


-4) 


1.7( 


-6) 


8.0( 


-5) 


3.6( 


-4) 


1.0( 


-2) 


2.4(-8) 




9. 3 (-9) 




<6.7(-10) 


Sco 


0.87 


1.4(21) 


1.4( 


-4) 


1.0( 


-7) 


6.5( 


-5) 


2.8( 


-4) 


1.0( 


-1) 


4.2(-9) 


3.1(-9) 


1.5(-9) 


<9.2(-10) 


<1.0(-10) 



^ Fractional abundance of species Y. 
Ntotai= A^(H I) + 2iV(H2). 
2.7(20)= 2.7x10^0. 



Table 8. Chemical Network 



Reaction Rate coefficient or Rate Notes 



1. 
2. 

3. 
4. 
5. 
6. 

7. 

8. 
9. 

10. 

11. 

12. 

13. 
14. 
15. 
16. 

17. 
18. 
19. 
20. 
21. 
22. 
23. 
24. 
25. 

26. 

27. 

28. 
29. 

30. 



C0+ + H ^ CO + H+ 
HCO+ + e- ^ CO + H 

CO + He+ ^ C+ + O + He 

CO + hi/-^C + 
CH+ + O ^ C0+ + H 
C0+ + H2 ^ HCO+ + H 

C0+ + e- ^ C + O 

H2 + He+ ^ He + H+ 
He + CR ^ He+ + e" 

He+ + e~ ^ He + h 1/ 

CH+ + e- ^ CH + H 

CH+ + e- ^ CH + H2 
C + H^CH + h:/ 
C+ + CH ^ CH+ + C 
C+ + CH ^ C^ + H 
C + CH ^ C2 + H 

CH + O ^ CO + H 
CH + h zy ^ CH+ + e" 
CH + hi^^C + H 
CH+ + H2 ^ CH+ + H 



H2 
O - 
c- - 
O - 
e- - 



^ CH+ 
HCO+ 



CH+ 

CH2+ 

cat 
CH3+ 

CH+ 
CH+ + e 

C+ + H2 

C+ + OH ^ C0+ + 
C+ + H2O ^ HCO+ 

CH + H ^ C + H2 



H 

hH 



H2 



C + H2 
HCO+ H 
CH2 +H 

C + H2 + H 
CHj + h 1/ 



H 



H 



7.5 X 10"^° cm^ 

/ o„„ \ 0-69 

2.4 X 10-7 X (^) cm3 s-i 
1.6 X 10-9 cm^ 

3.5 X 10-1° cm3 s-i 



1.4 X 10-9 cm^ 
1.0 X 10-7 X (^)°'^' cm3 s-i 
1.0 X 10-1^ cm^ 
1.1/n X 10-17 cm^ s-i 

/ N 0.64 

2.4 X 10-12 X (^j 
1.25 X 10-7 X 



cm"' s 
cm'' s 



2.45 X 10 



300 



cm s 



3 „-i 



1.0 X 10-17 cm^ s-i 
1.5 X 10-9 cm3 s-i 
5.4 X 10-1° cm^ s-i 
1.0 X 10-12 X s/Teff cm^ s-i 

9.5 X 10-11 X ^ 
4.2 X IQ-i" 



— ^^^^ cm s 
300 ^^^^ ^ 



X c 



-Tuv g-l 



1.3 X 10^9 X e-'"''^ s-i 
1.2 X 10-9 cm^ s-i 
1.2 X 10-9 cm^ s-i 
7.5 X 10-1° cm3 s-i 



1.25 X 10- 



^ cm3 s-i 
4.40 X 10^1" cm^ s-i 



7.0 X 10-^ X cm^ s-i 



3.5 X 10-s X cm^ 3-1 

cm^ s-i 

7.7 X 10-1° cm3 s-i 
9.0 X 10~i° cm^ s-i 



1.0 X 10-1^ 



1 
2 
1 

3 
3 

2 

3 
1 

4 

2 

2 

4 

5,6 
5,6 
1 



5.0 X 10-11 X 



3^ cm3 s 1 



Notes. - (1) Draine & Katz 1986; (2) Mitchell 1990; (3) Anicich 1993; (4) Millar et al. 
1997; (5) Fcdcrman & Huntress 1989; (6) Pan et al. 2001; (7) Barsuhn & Nesbit 1978; (8) 

Fedcrman ct al. 1997. 

^ G(CO) in units of 10-i° s"! are 1.6, 1.3, 1.2, 1.1, 1.1, 1.2, 1.4, 1.3, 1.1, 0.89, 1.6, and 
1.1 for T] Tau, e Per, A Ori, HD 112244 A, HD 112244 B, 1 Sco, tt Sco, 5 Sco, lo^ Sco, a Sco, 
23 Ori, and Sco, respectively. 



-25- 



Table 9. CH Fractional Abundances 



Star 


Observed" 


Model 1'' 




Model 2'= 




Model 3'' 




r] Tau 




<6. 1x10-12 


<2.8xl0-9 




<2. 7x10-9 




e Per 


<6. 2x10-9 


(2.2-5.3) X 10" 


ll 


(0.5-1.2) X 10- 


-8 


(4.2-7.4) X 10- 


-9 


A Ori 




(1.2-1.8)xl0- 


11 


(0.9-1.4) x 10- 


-9 


(4.6-5.2) x 10- 


10 


HD 112244 A 


4.0x10-'^ 


(3.3-5.4) X 10- 


11 


(2.2-3.6) X 10- 


-8 


(1.6-2. 2)xl0- 


-8 


HD 112244 B 


2.8x10-9 


(3.6-6.9) X 10- 


11 


(0.7-1.4) X 10- 


-8 


(5.2-7.4) x 10- 


-9 


1 Sco 


<1. 1x10-9 


(6.5-9.9) X 10" 


12 


(3.5-5.3) X 10- 


10 


(1.8-2.0)xl0- 


10 


TT Sco 


<8.9xl0-i° 


(0.6-1.3) X 10- 


11 


(0.8-1.8) X 10- 


-9 


(0.7-1. l)x 10- 


-9 


S Sco 


1.6x10-9 


(1.4-2.3)xl0- 


11 


(1.9-3. 0)xl0- 


-9 


(1.0-1. l)xl0- 


-9 


ui^ Sco 


1.9x10-9 


(3.7-5.4) X 10- 


11 


(0.9-1.3) X 10- 


-8 


(5.3-6.1)xl0- 


-9 


a Sco 


1.3x10-9 


(1.9-2.5)xl0- 


11 


(3.2-4.2) X 10- 


-9 


(1.5-1.6)xl0- 


-9 


23 Ori 


9.3x10-9 


(1.2-2.9)xl0- 


12 


(1.0-2.4) X 10- 


-9 


(0.8-1.5) x 10- 


-9 


Sco 


1.5x10-9 


(2.4-4.2) X lo- 


11 


(0.7-1.2) X 10- 


-8 


(4.2-5.6) x 10- 


-9 



^ References - Federman 1982; Danks ct al. 1984; Federman ot al. 1997. 
^ Equilibrium chemistry {vion = 0, Vneutrai = 0, a;(CH-'-) = 3.3x10^1^). 
^ Non-equilibrium chemistry {vion = 3 km s-"^, Vneutral = 0, observed a:;(CH+)). 
Non-equilibrium chemistry {vion = 3 km s-^, Vneutral = 3 km s-^, observed a;(CH-'-)). 
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Table 10. CO Fractional Abundances 



Star 


Observed" 


Model I'' 




Model 2" 




Model 3'' 




77 Tau 


<8.1xl0-^ 


<7.7xl0-i3 


<1. 6x10-10 


<1. 7x10-1° 


e Per 


2.7x10-9 


(2.8-6.9) X 10- 


12 


(2.4-4.6) X 10" 


10 


(2.6-5.8) X 10- 


10 


A Ori 


6.2x10-9 


(1.2-3.2) X 10" 


-11 


(1.6-4.0) X 10- 


10 


(1.8-4.4) X 10- 


10 


HD 112244 A 


1.8x10-8 


(4.0-7.6) X 10- 


-12 


(l.l-1.8)xl0- 


-9 


(1.2-2.3)xl0- 


-9 


HD 112244 B 


5.1x10-9 


(0.5-1.3) X 10- 


-11 


(3.8-9.8) X 10- 


10 


(0.4-1. 2)xl0- 


-9 


1 Sco 


<1.6xl0-9 


(0.4-1.7) X 10- 


-10 


(0.3-1.2) X 10" 


-9 


(0.4-1.3)xl0- 


-9 


TT Sco 


1.5x10-9 


(2.4-5.7) X 10- 


-12 


(0.8-1.6) X 10- 


10 


(0.8-1.8) X 10- 


10 


S Sco 


2.2x10-9 


(2.1-5.9)xl0- 


-11 


(0.5-1.3)xl0- 


-9 


(0.6-1.5)xl0- 


-9 


co-^ Sco 


5.3x10-9 


(l.l-2.1)xl0- 


-11 


(1.0-1. 7)xl0- 


-9 


(l.l-2.0)xl0- 


-9 


a Sco 


1.3x10-9 


(3.8-9.3) X 10- 


-11 


(1.5-3.4) X 10- 


-9 


(1.7-3.9)xl0- 


-9 


23 Ori 




(2.0-5.9)xl0- 


-12 


(1.7-4.9)xl0- 


10 


(1.7-5.0)xl0- 


10 


Sco 


3.1x10-9 


(0.8-4.7) X 10- 


-11 


(0.5-2.9) X 10- 


-9 


(0.6-3.4) X 10- 


-9 



^ Crcnny & Fcdcrman, private communication. 
Equilibrium chemistry {vion = 0, Vneutrai = 

0, a;(CH+) = 3.3x10-12). 
^ Non-equilibrium chemistry {vion = 3 km s-^, Vneutrai = 0, observed x(CH+)). 
Non-equilibrium chemistry {vion = 3 km s-^, Vneutrai = 3 km s~^, observed a;(CH-'-)). 
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Table 11. Wx (mA) of C I Lines toward rj Tau 



Wavelength 


Measured 


Calculated 


(A) 






J = 0: 






1260.736 


8.39 ± 1.51 


10.64 


1277.245 


16.38 ± 1.24 


19.88 


1280.135 


7.30 ± 1.09 


7.38 


1328.833 


16.59 ± 0.99 


15.46 


1560.309 


24.72 ± 0.86 


24.50 


1656.928 


33.87 ± 1.06 


34.47 
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Table 12. Wx (mA) of C 1 Lines toward e Per 



A Measured Previous Measurements Calculated 

(A) HST Copernicus JS« JJiT 



J = 0: 



1260.736 


18.77 ± 0.34 






14.0 ± 0.7 


17.32 


1328.833 




22.3 ± 0.6 




29.0 ± 1.5 


24.06 


J = 1: 












1260.927 


< 0.32 






< 1.3 


1.14 


1260.996 


< 0.32 






< 1-2 


0.86 


1261.122 


< 0.32 






< 1.2 


1.42 


1329.085 










1.84 


1329.100 




5.5^ ± 0.6 


3.5^ ± 1.8 


5.4^ ± 0.9 


2.28 


1329.123 










1.40 


J = 2: 












1261.426 


< 0.32 






< 1.4 




1261.552 


< 0.32 






< 1.4 




1329.578 




< 0.6'^ 


< 4.0<^ 


< 2.2'' 





^ Jenkins & Shaya 1979. 

^ Jenkins et al. 1983. 
Total Wx for AA1329.085, 1329.100, 1329.123. 

Total Wx for AA1329.578, 1329.601. 
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Table 13. Wx (mA) of C I Lines toward A Ori 



\ 

A 


Measured 


Other Measurements 






HR» 


MR« 




JJL^ 




3 = 0: 
1157.910A'' 


26.97 ± 0.26 








28.09 


1157.910B'' 


2.59 ± 0.30 








4.25 


1158.324A 


15.28 ± 0.23 








15.52 


1158.324B 


2.99 ± 0.35 








1.52 


1260.736A 




35.44 ± 0.33 


37.0^ ± 2.0 


38.0^ ± 2.6 


31.27 


1260.736B 




3.17 ± 0.24 






4.88 


1328.833 




48.23 ± 0.32 


52.0*= ± 4.0 


51.0*= ± 2.6 


45.66 


J = 1: 












1157.770 


5.82 ± 0.21 








5.58 


1158.035 


8.95 ± 0.21 








9.56 


1158.130 


4.42^ ± 0.21 








4.28 


1158.674 


0.85 ± 0.09 








1.25 


1158.732 


1.39 ± 0.10 








1.47 


1260.927 




5.17 ± 0.22 




< 6.3 


8.71 


1260.996 




9.00 ± O.oD 




D.O ± 2.1 


6.84 


1261.122 




9.73 ± 0.30 


13.0 ± 1.0 


7.6 ± 2.1 


10.43 


1329.085 










12.96 


1329.100 




47.839 ± 0.32 


42.09 ± 2.0 


56.09 ± 3.0 


15.25 


1329.123 










10.39 


J = 2: 












1158.019 


2.15 ± 0.21 








4.03 


1158.132 










0.44 


1158.397 


0.87 ± 0.21 








0.74 


1261.426 




1.95 ± 0.29 


3.4'^ ± 1.2 


< 6.3 


1.45 


1261.552 




4.82 ± 0.32 




< 6.6 


4.17 


1329.578 




11.68* ± 0.32 


14.0* ± 2.0 


9.5' ± 3.7 


6.60 


1329.601 










2.34 



* HR: High Resolution spectra, MR: Medium Resolution spectra. 

Jenkins & Shaya 1979. 
^ Jenkins et al. 1983. 

^ A: Main component B: Secondary component, if resolved. 
^ Involves both the main and secondary J = components, 
f Blended with A1158.132 (J = 2). 
s Total Wa for AA1329.085, 1329.100, 1329.123. 

1^ Total Wx for AA1261.426, 1261.552. 
' Total Wx for AA1329.578, 1329.601. 
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Table 14. (mA) of C I Lines toward HD 112244 



A 


Measured 


Calculated 


Component A 


Component B 


Component A 


Component B 


J = 0: 
1188.833" 


53.66 ± 2.42 


32.44 ± 2.05 


29.26 


25.16 


1193.031'' 


52.65 ± 2.46 


47.22 ± 2.51 


45.28 


39.42 


1193.996 


22.36 ± 1.91 


25.10 ± 2.24 


24.77 


21.23 


1260.736 


39.11 ± 1.20 


27.60 ± 1.26 


48.07 


41.85 


1276.483 


10.16 ± 0.68 


14.27 ± 0.98 


12.31 


10.46 


1277.245^= 


75.27 ± 1.28 


63.47 ± 1.16 


61.93 


54.31 


1280.135 


44.19 ± 1.40 


29.04 ± 1.07 


40.22 


34.80 


1560.309 


76.95 ± 2.71 


61.83 ± 2.35 


75.89 


66.56 


J = 1: 










1188.993" 


13.37 ± 


1.77 


2.13 


2.12 


1193.009'' 






9.54 


9.90 


1193.264 






4.06 


4.08 


1193.679 






3.72 


3.73 


1260.927 






5.22 


5.26 


1260.996 






4.04 


4.05 


1261.122 






6.34 


6.44 


1277.282'= 






18.42 


20.23 


1277.513 


28.50 ± 


0.99 


8.92 


9.19 


1279.890 


0.40 ± i.io 




5.57 


5. DO 


1560.682'^ 


32.71 ± 2.35 


35.82 ± 2.26 


22.69 


24.95 


1560.709'' 






11.03 


11.37 


J = 2: 










1193.240 






1.96 


2.81 


1193.393 






0.36 


0.53 


1193.649 






0.02 


0.04 


1261.426 






0.60 


0.87 


1261.552 






1.75 


2.52 


1277.550 






4.57 


6.38 


1277.723 






0.90 


1.31 


1561.340*= 






1.12 


1.62 


15G1.438'' 


17. GO ± 


2.21 


r).()7 


7.91 



^ Neutral carbon lines are obscured by a strong CI I. 

^ Blended lines. The total W\ is listed for the stronger line. 

*= Blended lines. The total Wx is listed for the stronger line. 

Blended lines. The total W\ is listed for the stronger line. 
*= Blended lines. The total W\ is listed for the stronger line. 



Table 15. Wx (mA) of C I Lines toward 1 Sco 



A M(>;isurc(l Other iiicasur(>iiK>ut Calculated 

(A) EB^ MR^ JJiT 



J = 0: 
1193.031A'= 
1193.031B 
1193.996A<^ 
1193.996B'' 
1260.736 
1277.245 
1280. 135A 
1280. 135B 
1328.833A 
1328.833B 
1560.309A 
1560.309B 
1656.928A 
1656.928B 

J = 1: 
1193.009^ 
1193.679 
1194.406 
1260.927 
1260.996 
1261.122 
1277.5139 
1279.890 
1280.597 
1329.085 
1329.100 
1329.123 
1560.682 
1560.709 
1656.267 
1657.379 
1657.907 

J = 2: 
1261.426 
1261.552 



10.43 ± 0.62 
11.46 ± 0.62 
3.63 ± 0.42 
3.02 ± 0.62 



12.97 ± 0.66 
15.82 ± 0.98 



2.89 ± 0.46 
4.19 ± 0.52 
1.46 ± 0.37 



6.59 ± 0.52 
9.26 ± 0.56 
7.05 ± 0.63 



19.73^= ± 0.98 
69.20^ ± 1.42 
7.38 ± 0.70 
7.06 ± 0.73 
14.50 ± 0.86 
14.50 ± 0.85 
26.72 ± 0.83 
14.15 ± 0.75 
37.89 ± 0.97 
28.58 ± 0.90 



5.23 ± 0.94 
3.45 ± 0.78 
5.11 ± 0.81 
4.75 ± 1.15 
5.56 ± 0.74 
2.20 ± 0.90 

26.90'' ± 1.10 

32.4' ± 0.98 

25.89 ± 0.93 
19.09 ± 0.86 
22.75 ± 0.94 

1.20 ± 0.94 
3.68 ± 0.68 



22.0^= 



34.0*= 



± 2.4 



± 4.4 



< 6 

< 6, 

< 6, 



38.0^ 



± 3.3 



< 6 

< 6, 



9.35 
9.56 
3.62 

3.37 
20.18*= 
37.48*= 
7.11 
6.94 
13.98 
15.23 
21.13 
25.06 
28.11 
37.20 

10.35 
4.21 
1.67 
5.85 
4.57 
7.05 
9.76 
6.23 
3.39 
8.82 
10.47 
7.01 
23.87 
12.07 
25.70 
19.42 
22.96 

2.11 
5.41 
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Table 15 — Continued 



A 

(A) 


Measured 


Other measurement 


Calculated 




MR" 


JJL" 


1277.550S 




13.07 ± 0.87 




10.89 


1277.723 




1.80 ± 0.89 




3.05 


1280.333 




2.50 ± 0.90 




3.09 


1329.578 


6.76 ± 0.55 


10.30-' ± 0.94 


< 8.7^' 


7.91 


1329.601 


4.63 ± 0.72 






3.30 


1561.340 




3.26 ± 0.67 




3.78 


1561.438 




13.08 ± 0.81 




13.42 


1657.008 




15.43 ± 0.75 




18.04 


1658.121 




8.99 ± 0.68 




9.91 



HR: High Resolution spectra, MR: Medium Resolution spectra. 
^ Jenkins et al. 1983. 
Blended Hnes. 

A: Main component B: Secondary component, if resolved. 
° Involves both the main and secondary J = components. 
^ Involves a second component and A1277.282. 
s Blended lines. 

^ Total Wx for AA1329.085, 1329.100, 1329.123. 
' Total Wx for AA1560.682, 1560.709. 
j Total Wx for AA1329.578, 1329.601. 
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Table 16. Wx (mA) of C I Lines toward tt Sco 





A/T API C!1 1 Ti^ri 


(l^" rn^f A/Tap qi 1 fi^TTi i" 




(A.) 








J = 0: 








1260.736 


3.40 ± 0.95 


4.9 ± 2.0 


4.84 


1277.245'' 


13.46 ± 0.49 




11.28 


1280.135 


2.80 ± 0.50 




3.14 


1328.833 


8.97 ± 0.32 




7.69 


1560.309 


10.63 ± 0.66 




13.97 


1656.928 


26.74 ± 0.74 




24.74 


J = 1: 








1260.927 


< 0.94 




0.25 


1260.996 


< 0.94 




0.19 


1261.122 


< 0.94 


< 0.7 


0.32 


1277.282 






1.40 


1277.513 


< 0.38 




0.48 


1279.890 


< 0.50 




0.27 


1280.404 


< 0.50 




0.09 


1329.085 






0.41 


1329.100 


1.80^ ± 0.32 




0.52 


1329.123 






0.31 


1560.682 


2.17'' ± 0.64 




1.73 


1560.709 






0.59 


1656.267 


1.92 ± 0.50 




1.89 


1657.379 






1.15 


1657.907 






1.53 


J = 2: 








1261.552 


< 0.94 


< 0.8 





Jenkins & Shaya 1979. 
^ Involves A1277.282. 

^ Total Wx for AA1329.085, 1329.100, 1329.123. 
^ Total Wx for AA1560.682, 1560.709. 



-34- 



Table 17. W\ (mA) of C I Lines toward 6 Sco 



A 

(A) 


Measured 


Other measurement 


Calculated 


TTT> O. 


MR" 




J = 0: 

1157. giOA'^ 


23 06 ± 28 






20.36 


1157.910B'= 


2 11 ± 20 






4.87 


1158.324A 


10 84 ± 23 






10.78 


1158.324B 


1.44 ± 0.25 






1.89 


1260. 736A 


21.46 ± 0.43 


28.53"^ ± 0.22 


32.0"^ ± 3.0 


22.72 


1260. 736B 


6 03 ± 36 






5.56 


J = 1: 










1157.770 


5.06 ± 0.21 






5.38 


1 1 58 0S5'° 


7.98 ± 0.21 






8.51 


1158.130-^ 


4.05 ± 0.22 






4.24 


1158.674 


1.37 ± 0.14 






1.32 


1158.732 


1.59 ± 0.13 






1.53 


1260.927 


7.63 ± 0.30 


8.00 ± 0.20 




7.99 


1260.996 


6.54 ± 0.31 


6.50 ± 0.19 


7.2 ± 1.0 


6.49 


1261.122 


8.15 ± 0.31 


8.80 ± 0.21 


9.7 ± 2.0 


9.27 


J ^2: 










1158.019^ 


4.06 ± 0.26 






3.94 


1158.132^ 


0.43 ± 0.22 






0.49 


1158.397 


0.81 ± 0.21 






0.81 


1158.967 


0.45 ± 0.09 








1261.426 


1.61 ± 0.30 


2.74 ± 0.26 




1.55 


1261.552 


3.87 ± 0.35 


4.95 ± 0.21 


4.6 ± 1.4 


4.09 



HR: High Resolution spectra, MR: Medium Resolution spectra. 
^ Jenkins & Shaya 1979. 

A: Main component B: Secondary component, if resolved. 
^ Involves both the main and secondary J = components. 
® Blended lines. 
^ Blended lines. 
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Table 18. Wx (mA) of C I Lines toward Sco 



A 

(A) 


Measured 


Other measurements 


Calculated 


HbT Copernicus 




JJL" 


J = 0: 










1260.736 


35.54 ±0.14 


26.0 ± 1.7 


33.0 ± 3.6 


33.16 


1276.483 


6.6 ±4.1 




< 6.4 


8.85 


1328.833 


22.2 ±4.1 




62.0 ± 6.2 


39.63 


J = 1: 










1260.927 


8.19 ± 0.12 


4.7 ± 1.6 


< 5.8 


6.46 


1260.996 


6.42 ± 0.11 


2.7 ± 1.6 


< 5.7 


5.31 


1261.122 


8.22 ± 0.13 


11.0 ± 1.7 


< 5.7 


7.41 


1329.085 








8.80 


1329.100 


42.5^ ± 4.1 




47.0^ ± 5.2 


9.81 


1329.123 








7.49 


J = 2: 










1261.426 




< 3.3 


< 5.8 


1.31 


1261.552 


4.34 ± 0.14 


< 3.3 


< 6.6 


3.46 


1329.578 


4.9"^ ±4.1 




< 7.1** 


5.16 


1329.601 








2.06 



Snow & Jenkins 1980. 
^ Jenkins et al. 1983. 

" Total Wx for AA1329.085, 1329.100, 1329.123. 
Total Wx for AA1329.578, 1329.601. 
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Table 19. Wx (mA) of C I Lines toward a Sco 



A 

(A) 


Measured 


Other measurement 


Calculated 


HR" 






T (\. 

J — U. 


23.59 ± 0.64 


36.08"^ ± 0.50 






IZUU. ( OU-D 


2.41 ± 0.41 






A Ar\ 


±0^0.000 




50.90'' ± 0.50 




Oiy.O i 


xouu.out? 




56.50^^ ± 0.54 




Q7 


J = 1: 










1260.927 


12.46 ± 0.51 


9.40 ± 0.50 




11.90 


1260.996 


10.04 ± 0.47 


11.10 ± 0.50 




9.83 


1261.122 


11.87 ± 0.48 


11.84 ± 0.50 


14.0 ± 2.0 


13.61 


1329.085 








16.12 


1329.100 




42.80* ± 0.50 




17.91 


1329.123 








13.78 


1560.682 




28.60 ± 0.45 




29.34 


1560.709 




19.10 ± 0.38 




20.80 


J = 2: 
1261.426 




3.60 ± 0.50 


4.3 ± 2.0 


3.11 


1261.552 




6.50 ± 0.50 


11.5 ± 2.0 


7.53 


1329.578 




9.40 ± 0.37 




10.55 


1329.601 




6.20 ± 0.37 




4.78 


1561.340 




1.69 ± 0.32 




5.50 


1561.438 




16.10 ± 0.43 




16.69 



HR: High Resolution spectra, MR: Medium Resolution spectra. 
^ Jenkins & Shaya 1979. 

A: Main component B: Secondary component, if resolved. 

Involves both the main and secondary J = components. 
" Total Wx for AA1329.085, 1329.100, 1329.123. 
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Fig. 1. — Observed (histograms) and model (lines) spectra of the C I multiplets toward 
1 Sco. All spectra are medium resolution unless otherwise noted, a) The top, middle, and 
bottom rows from left to right show multiplets AA1193, 1194 (high resolution), 1260, 1277, 
1280, and 1329 (high resolution, J — and 1 levels only), respectively. The residuals are 
displayed above each spectra, offset by +1.15. 
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Fig. 1. — b) Spectra of multiplets AA1329 (high resolution, J — 2 level only), 1329, 1561, 
and 1657. 



-43- 



200 
150 
100 
50 


I 

200 



, 

\ 

\ 


1 

ri Tau 


\ 

\ 

\ 

\ 





200 



150 



100 




10 20 30 40 50 10 20 30 40 50 



4> 150 
U 

S 

CQ 100 

u 

s 

H 



\ 


lOri 

\ \ 
\ \ 






.... 1 .... 1 .... 1 ... . 



200 



100 



50 





\\ 


HD 112244 A 




\\ 


\\ 


. 1 . . , 1 . . . 1 . . . 







50 



200 
150 
100 
50 




100 150 200 

200 



150 



1 


1\ HD 112244 B 

\\\^\ ; 




\ \ ^ 












1 ... 1 . 





100 



50 



20 



40 



60 





80 



1 


i 

1 SCO 

\\ 




\ \^ 




. . . , 1 . . . , 1 . . . . 1 . . , 



400 800 



1200 



Density (cm'^) 

Fig. 2. — The density-temperature relationships for the main components along our sight 
lines. The curves for the secondary component toward HD 112244 (component B) are also 
displayed. The two solid lines represent the density-temperature correspondence calculated 
from the lowest and the highest nj=i/nj=o ratios allowed by our C I column densities. A sim- 
ilar pair of lines (dashed lines) show the relationships calculated from the nj=2/'nj=o ratios. 
If only an upper limit was available for either the J = 1 or 2 column density, then a single 
line shows the maximum gas densities as a function of temperature for the corresponding 
ratio. The density-temperature values, inferred from a ratio, are between the curves be- 
longing to the same ratio. In the cases of upper limits, all density-temperature points that 
are to the left of the curve are allowed. The dotted horizontal lines are the H2 excitation 
temperatures from Savage et al. 1977. The thick rectangles contain all density-temperature 
values that are consistent with both ratios and the H2 excitation temperatures. The overlap 
between the ranges was not sufficient toward HD 112244 B, 1 Sco, and Sco; in these cases, 
the rectangles contain the extreme density-temperature values that are consistent with all 
measures. 
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Fig. 2.— Cont. 
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Fig. 3. — CH produced by the non-equilibrium CH+ chemistry, normahzed by the observed 
amount of CH+, as a logarithmic function of CH"^ columns. The arrows and triangles indicate 
the upper limits and measured amount of CH, respectively. The maximum and the minimum 
amount of CH associated with the CH+ chemistry are shown by ranges. The identifier of 
each sight hne is printed above the symbols except for Sco which is below the symbol for 
clarity. 



